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ABSTRACT 

We study the previously discovered hemispherical power asymmetry (Eriksen et al. 2004, Hansen et 
al. 2004, Hansen et al. 2009) with the WMAP 9-year data. We find that the asymmetry is statistically 
significant at the 3.2<j confidence level for £ = 2-600, where the data is signal dominated, using the 
WMAP 9-year KQ85 mask. The preferred asymmetry direction is (Z,6) = (227,-27), 44° away 
from the Ecliptic South Pole. Individual asymmetry axes estimated from six independent multipole 
ranges are all consistent with this direction. This agreement across a wide range of scales argues 
against a statistical a-posterior 1 interpretation of the effect. We estimate cosmological parameters on 
different parts of the sky and show that the parameters A 3 , n s and are the most sensitive to the 
asymmetrically distributed power asymmetry. In particular, for the two opposite hemispheres aligned 
with the preferred asymmetry axis, we find n s = 0.96 ± 0.022 and n s = 0.99 ± 0.024, respectively. 
Subject headings: cosmic microwave background — cosmology: observations — methods: statistical 
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1. INTRODUCTION 

Sho rtly after the releas e of the first-y e ar W MAP 
data (iBennett et all 12003( 1. lEriksen et atl ([2004D and 
lHansen et all (|2004f ) reported a detection of a hemispher- 
ical power asymmetry in the cosmic microwave back- 
ground (CMB) on large angular scales, £ = 2-40. The 
power in this multipole range was found to be signif- 
icantly higher in the direction towards Galactic longi- 
tude and latitude (Z = 237°, b = -20°) than in the op- 
posite direction. Due to computational limitations at 
the time, higher multipoles were not investigated. These 
findin gs we r e supported by numerous other studies, e.g., 
iParkl (|2004f ): lHansen et all (|2009f ) and references therein. 
However, the significance of the results has often been 
called into question, in particular, due to the alleged a- 
posteriori nature of the statistics used. In particular, it is 
debated whether the statistic has been designed to focus 
on visually anom alous features reveal ed by an inspection 
of the data (e.g.. iBennett et al.l[20lTh . 

The only rigorous way to contend with this assertion 
is by performing repeated experiments and analysing the 
resulting independent data sets that may provide addi- 
tional information. For cosmological studies, this is in 
general difficult, given that there is only one available 
Universe. However, it is not impossible — the stan- 
dard inflationary cosmological model assumes that the 
Universe is homogeneous and isotropic, and that the ini- 
tial fluctuations have amplitudes that are Gaussian dis- 
tributed, independent and with random phase. This im- 
plies that different physical scales should be statistically 
uncorrelated, and therefore the morphology of the largest 
scales should not have any predictive power over the mor- 
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phology of the smaller scales. For the power asymmetry, 
this suggests that there is a possibility to study effec- 
tively new data sets by considering angular scales that 
have not previously been studied. 

This exten sion to smalle r angular scales was first un- 
dertaken by lHansen et al.l ([200 9). when analyzing the 
WMAP 5-year temperature data set. The asymmetry 
was then found to extend over the range £ = 2 — 600 
with a preferred direction (I = 226°, b = —17°) for the 
higher multipoles, fully consistent with the direction for 
the lower multipoles found in the original 1-year WMAP 
analysis. Two approaches were used for the analysis: (1) 
a statistical model selection procedure taking into ac- 
count the penalty for including 3 new parameters (am- 
plitude and direction of asymmetry), which showed that 
indeed an asymmetric model was preferred; and (2) a 
simple test in which the preferred power asymmetry axis 
was estimated independently for six multipole bins of 
width A£ = 100. It was found that these directions, 
which should be statistically independent, were strongly 
aligned; none of the 10 000 simulated isotropic CMB 
maps showed a similarly strong clustering of preferred 
directions. An alternative approach modeled the power 
asymmetry i n terms of a dipola r modulation fi e ld, as 
suggested bv lGordon et all (|20Q5h . iHoftuft et all ([20091 ) 
found a 3.3<j detection using data smoothed to an angu- 
lar resolution of 4.5° FWHM, with an axis in excellent 
agreement with previous results. These studies, covering 
very different angular scales than those used in the origi- 
nal analysis, argue against an a-posteriori' interpretation 
of the effect. 

In this Let t er, w e repeat the high-^ analysis due to 
lHansen et all ([2009D using the WMAP 9-year (hereafter 
WMAP9) data. However, the main goal is to estimate 
cosmological parameters in the two maximally asymmet- 
ric hemispheres, in order to assess their stability with 
respect to the power asym metry (for a close ly related 
theoretical study, see, e.g., Moss et all ([20111 ) and ref- 
erences there i n). A similar analysis was performed in 
lHansen et all (|2004f ) using the WMAP first-year data, 
but only taking into account the asymmetry observed in 
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Fig. 1. — Regions used for parameter estimation. 
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Fig. 2. — Dipole directions for six 100-multipole bins of the local 
power spectrum distribution derived from the WMAP9 combined 
V- and W-band map. We also show the direction for the full t = 
2 — 600 range (white cross), for the £ = 2 — 40 interval determined 
from the WMAP first year data (green circle) and the £ = 2 — 600 
range from WMAP5 (black cross). NEP and SEP denotes the 
North and South Ecliptic Poles. The dipole directions for the local 
parameter estimate maps are also shown. 

the t = 2 — 40 range, limited by a grid-based approach, 
and consequently only considering a few cosmological pa- 
rameters. In the following, we use CosmoMC to obtain 
the full posterior of all relevant cosmological parameters 
using the entire multipole range afforded by the WMAP9 
data. 

2. DATA AND METHOD 

We use the publicly avail able WMAP 9 tempe rature 
sky maps (Bennett et al. iBennett et al.l (|2Q12f )). co- 
adding (with inverse-noise-variance weighting) the V 
(61 GHz) and W (94 GHz) band foreground-cleaned 
maps. We also generate a set of 10 000 simulated CMB- 
plus-noise Monte Carlo (MC) simulation s based on the 
WMA P best-fit ACDM power spectrum Hin shaw et al.l 
([2012h . noise rms maps and beam profiles for the V and 
W bands. The WMAP9 KQ85 Galactic and point source 
mask is used to remove pixels with high foreground con- 
tamination. 

Power asymmetry: — We estimate the Ct auto power 
spectra locally on (1) regions corresponding to the twelve 
pixels of an = 1 HEALPix map (always applying 

the Galactic and point source mask), illustrated in Fig- 
ure [TJ (2) on the two opposite hemispheres defined by 



the direction of maximum asymmetry; (3) on disks of di- 
ameter 90 degrees centered o n the directions of maximum 
asymmetry. The MASTER (jHivon et al.ll2002D approach 
is used to estimate the power spectra. When applying 
the MASTER kernel to determine the pseudo-spectra, 
we bin the spectra into bins of width A£ = 16. We apply 
the same procedure to the set of 10 000 simulated maps. 

To estimate the dipole directions of the local spectra, 
we combine the 16^-bins further into blocks containing 
ap proximately 100 mult ipoles, following the binning used 
in lHansen et al.l (|2QQ9h . For each 100-multipole block 
determined at a given pixel, the mean power computed 
from the simulations is subtracted, and the results are 
divided by the similarly evaluated standard deviation in 
order to prevent the directions close to the Galaxy, where 
the Galactic mask increases the sample variance, domi- 
nating the fits due to possible large fluctuations. 

In this way, we obtain an N s ^ e = 1 map for each 
100-multipole block as described above, and the value in 
each pixel shows the normalized deviation away from the 
best fit full-sky spectrum. We make a spherical harmonic 
transform of this map, and extract the dipole amplitude 
and asymmetry direction. 

For an isotropic map, the power spectrum should be 
uncorrelated between multipoles. Although masking 
does introduce correlations between adjacent multipoles, 
it is not expected that this will be significant between the 
100-multipole blocks, and therefore the dipole directions 
should be random. The degree of alignment between the 
dipole directions of different multipole blocks is then used 
as a measure of power spectrum asymmetry. 

Cosmological parameter estimation: — Our primary inter- 
est here is to evaluate the directional dependence of the 
cosmological parameters in the temperature data. We 
apply a similar binning and power spectrum estimation 
method as described above. For the high-^ likelihood, 
we compute the MASTER estimates of the spectra from 
the full-resolution A" S id e = 512 map. For the \ow-£ likeli- 
hood, we replace the WMAP pixel likelihood by a MAS- 
TER estimate computed for a single bin, £ G [2, 31]. The 
maximum multipoles used in the parameter analysis are 
f max = 1008 for the analysis on regions defined by the 
Aside = 1 pixels, and £ max = 608 for the hemispheres 
defined by the direction of maximum power asymmetry. 
For all spectra, w e subtract the best-fi t unresolved point 
source amplitude (Hinshaw et al. 2012) before parameter 
estimation. 

Our TT likelihood code uses the offset log-normal 
term that was intr oduced into the WMAP likelihood in 
I Verde et~aLl (|2003l ) - hence the total likelihood is a linear 
combination of Gaussian and log- normal terms: 

- logC(C b \C b ) ~ 1 £ AC b C b - b } AC£+| Az bQ bb '^4 

b,b' b,b' 

where Cb and Cb are the estimated and model power 
spectra respectively, ACb = C5 — C5, C bb } is the covariance 
matrix, estimated using CMB plus noise MC simulations, 
Zb = ln(C6 + A5), (where A5 is the noise spectrum) and 

Qbb' = (Cb + Nb)C^ b } (Cb' + Nv) is the local transforma- 
tion of the covariance matrix to the log-normal variables 
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Fig. 3. — Upper: Power spectra estimated from two hemispheres defined by the preferred dipole direction determined from the £ - 
range, as computed within disks of diameter 90°. Lower: Difference between the two power spectra. 
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Fig. 4. — Estimated ACDM parameters: in red the WMAP9 full sky estimates, and in black estimates within regions defined by an 
^side — 1 map. The values on the x-axes correspond to pixel number on the map. The estimates are based on multipoles in the range 
£ = 2 - 1008. 



To construct the covariance matrix, we use 500 CMB 
plus noise Gaussian simulations. The covariance matrix 
propagates the uncertainties introduced by effects such 
as the noise, mask geometry, and associated sample vari- 
ance. We make no attempt to include beam uncertainty 



in our pipeline. 

The fractional areas of the 12 A^de = 1 patches range 
from / s ky = 0.019 in the Galactic center to / s k y = 0.085 
in regions at high latitude. With such small patches, one 
might be concerned about the correctness of our likeli- 
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Fig. 5. — Top: Summary of n s -log(10 10 A s ) posterior in terms of 
la contours for the 12 regions. Middle: The same, but evaluated 
from two opposite hemispheres aligned with the preferred asym- 
metry direction. Colored dots indicate the values of Ho. Bottom: 
Same as the middle plot, but for Q h h 2 -\og(10 10 A s ) . 

hood approximation. We have confirmed that the pa- 
rameter estimates are unbiased, performing parameter 
estimation on some of the small regions in 500 simulated 
maps with known input parameters. 

The final posterior distribution is comprised of the 
product of the likelihood and an eventual prior distribu- 



tion that describes our previous knowledge of the param- 
eters. Since we do not consider polarization in our anal- 
ysis, we adopt the WMAP9 estimate of z re i = 10.6 ±1.1 
as a Gaussian prior on the re- ionization redshift. 

3. RESULTS 

In Figure [2j we show the first six dipole directions 
of the local power spectra as estimated on 12 pixels in 
100-multipole bands, together with the dipole for the 
total multipole range £ = 2 — 600. These directions 
are consistent with t hose found from the WMAP 1-year 
(Hansen et al.l 120041 ) and WMAP 5-year (jHansen et al.l 
|2009[ ) data sets, which are also indicated in the figure. 

For a statistically isotropic CMB temperature distribu- 
tion, the dipole directions from uncorrelated power spec- 
trum estimates should be distributed randomly on the 
sky. To quantify the significance of the power asymme- 
try, we consider the dispersion angle, which is defined as 
the mean angle, # meari , between all possible combinations 
of 100-multipole dipole directions up to a given ^ max - 
The expected dispersion angle for Gaussian simulations 
is 90 degrees, as confirmed by simulations. We calculate 
#mean(^max = 600) for the WMAP9 data and compare it 
to the distribution obtained from 10 000 CMB plus noise 
simulations. We found that only 27 of these exhibited a 
lower dispersion angle, implying a 3.2a significance for 
the powe r asymmetry. This is lower than previously re- 
ported in lHansen et al.l ([2009) , where none of the 10 000 
simulations had a similarly large mean angle. However, 
there are several changes in this analysis; (1) the new 
9-year Galactic and point source masks remove a larger 
fraction of the sky, thus increasing the scatter on the 
dipole directions due to increased sample variance; (2) 
larger disks are now used since the smaller disks from the 
5-year analysis result in several patches near the Galactic 
centre with an extremely small sky fraction in combina- 
tion with the new Galactic mask; and (3) 12 independent 
patches are now used instead of 3072 overlapping ones to 
speed-up the computations. It is also interesting to note 
that we have tried a number of permutations of the indi- 
vidual yearly sky maps and found variations depending 
on how many years are included, and which particular 
years are excluded. The results always remain highly 
significant, and the variations are found to be within the 
expected noise fluctuations as tested with simulations. 

In Figure [3j we show the difference between the power 
spectra computed for disks of diameter 90°, the so-called 
positive one centered at (Z,b) = (226°,— 27°), i.e., the 
dipole direction for the full t = 2-600 range, and its 
opposite (negative) disk. This plot visually illustrates 
the power asymmetry seen in the WMAP9 sky. Next we 
consider if this asymmetry in power is reflected in fits to 
the standard ACDM cosmological parameters. 

Figure [H shows the directional dependence of the six 
main ACDM parameters. The computed values and their 
standard deviations are shown in black, and the corre- 
sponding results from the WMAP9 analysis on the KQ85 
sky coverage is shown in red. The numbers on the x-axis 
correspond to the A^de = 1 pixel numbers. Inspecting 
the plots carefully, one finds that the majority of param- 
eter estimates fall within ~lcr of the WMAP9 full sky 
value. In Figure El we also show the dipole directions 
of the A^ide = 1 parameter maps. Clearly n s , A s and 
Qb seem to show a directional dependence similar to the 



5 



power spectrum asymmetry and these seem to be the pa- 
rameters mostly affected by the asymmetry. In the left 
top panel of Figure [5j we demonstrate the A s — n s corre- 
lation with la contours for each region. All contours are 
consistent with each other at better than 2<j, but some 
directional dependence is visible. 

We also estimated parameters using hemispheres cor- 
responding to the preferred power asymmetry direction 
for £ = 2 — 600. The first case we considered was re- 
stricted to £ max = 608 to cover only that part of the 
spectrum which is highly signal dominated and where 
the asymmetry is prominent. The absence of high mul- 
tipoles leads to large uncertainties in the parameters of 
interest and the asymmetry becomes less visible. We 
then repeat the analysis but with f max = 1008. In this 
case the error bars on the parameters are improved and 
the error ellipses for A s vs Vt^h 2 and A s vs n s computed 
on the positive (power-enhanced) and negative (power- 
deficit) hemispheres show a slight shift, as shown in the 
bottom panels of Figure [5] The best-fit parameters for 
each hemisphere lie just at the border or the la con- 
tours from the opposite hemisphere. The two maximally 
asymmetric hemispheres do not, therefore, indicate pa- 
rameter values significantly different from the WMAP9 
full-sky results. It is interesting to note that the power- 
deficit hemisphere prefers in general a higher Hq. The 
marginalised value obtained for the scalar spectral in- 
dex in the two hemispheres is n s = 0.959 ± 0.0224 and 
n s = 0.989 ± 0.0238 respectively (an w 1.3a difference). 
It is interesting to note that in one hemisphere, the spec- 
tral index is different from 1 at almost 2a whereas in the 
other it is fully consistent with 1. 

4. CONCLUSIONS 
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We measure a statistically significant power spectrum 
asymmetry in the WMAP9 temperature sky maps with 
3.2a significance as measured by the mean dispersion 
among the preferred directions derived from six (nearly) 
independent multipole ranges between t = 2 and 600, 
using the conservative mask adopted in the WMAP 9- 
year analysis. Only 27 out of 10 000 simulations show 
a similarly strong alignment. Conversely, the cosmo- 
logical parameters do not show a strong asymmetry, al- 
though the parameters A 3 , n s and ^5 do indicate some 
sensitivity to the hemispherical power asymmetry. The 
average preferred direction points toward Galactic co- 
ordinates (Z,6) = (226,-27), or Ecliptic coordinates of 
(Z,6) = (81, -46), i.e., 44° away from the South Ecliptic 
Pole. Thus, contrary to a common misconception, the 
CMB power asymmetry axis is not strongly aligned with 
the Ecliptic plane. 
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